THE CAUCHY PROBLEM FOR A FORCED HARMONIC OSCILLATOR 
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Abstract. We construct an explicit solution of the Cauchy initial value problem for the one- 
dimensional Schrodinger equation with a time-dependent Hamiltonian operator for the forced har- 
monic oscillator. The corresponding Green function (propagator) is derived with the help of the 
generalized Fourier transform and a relation with representations of the Heisenberg-Weyl group 
N (3) in a certain special case first, and then is extended to the general case. A three parameter 
extension of the classical Fourier integral is discussed as a by-product. Motion of a particle with a 
spin in uniform perpendicular magnetic and electric fields is considered as an application; a tran- 
sition amplitude between Landau levels is evaluated in terms of Charlier polynomials. In addition, 
we also solve an initial value problem to a similar diffusion-type equation. 



1. Introduction 

The time-dependent Schrodinger equation for the one-dimensional harmonic oscillator has the 
form 

th^ = Hi^, (1.1) 



where the Hamiltonian is 



2 \ dx"^ J 2 
Here a) and a are the creation and annihilation operators, respectively, given by 



H= — (-—+x')= — {aa^ + a^a). (1.2) 



. 1 f d\ 1 f d , , , 

a' = —;= [ X — — \ , a = —;= I X + — ; (1.3) 



a/2 V dx J ' a/2 V 9x ^ 

see [21] for another definition. They satisfy the familiar commutation relation 

[a, a^] = aa^ — a^a = 1. (1-4) 
A natural modification of the Hamiltonian operator (11. 2p is as follows 

H ^H{t) = ^ [aa) + a^a) +n{6{t)a + 6* (t) a^) , (1.5) 

where 5 [t) is a complex valued function of time t and the symbol * denotes complex conjugation. 
This operator is Hermitian, namely, (t) = H {t) . It corresponds to the case of the forced harmonic 
oscillator which is of interest in many advanced problems. Examples include polyatomic molecules 
in varying external fields, crystals through which an electron is passing and exciting the oscillator 



Date: February 1, 2008. 

1991 Mathematics Subject Classification. Primary 81Q05, 33D45, 35C05, 42A38; Secondary 81Q15, 20C35. 

Key words and phrases. The Cauchy initial value problem, the Schrodinger equation, forced harmonic oscillator. 
Landau levels, the hypergeometric functions, the Hermite polynomials, the Charlier polynomials. Green functions, 
Fourier transform and its generahzations, the Heisenberg-Weyl group N (3). 

1 



2 



RAQUEL M. LOPEZ AND SERGEI K. SUSLOV 



modes, and other interactions of the modes with external fields. It has particular applications in 
quantum electrodynamics because the electromagnetic field can be represented as a set of forced 
harmonic oscillators [9], [20], [23], [M], [35] and [15]. Extensively used propagator techniques were 
originally introduced by Richard Feynman in [16], [IT], [18] and [19] . 

On this note we construct an exact solution of the time-dependent Schrodinger equation 



where ipo (x) is an arbitrary square integrable complex valued function from {—oo, oo) . We shall 
start with a particular choice of the time-dependent function 6 {t) given by (13.61) below, which is 
later extended to the general case. The explicit form of equation fll.6p is given by (14. ip and (llU.ip 
below, and an extension to similar diffusion-type equations is also discussed. 

This paper is organized as follows. In section 2 we remind the reader about the textbook solution 
of the stationary Schrodinger equation for the one-dimentional simple harmonic oscillator. In sec- 
tion 3 we consider the eigenfunction expansion for the time-dependent Schrodinger equation (11.60 
and find its particular solutions in terms of the Charlier polynomials for certain forced harmonic 
oscillator. The series solution of the corresponding initial value problem is obtained in section 4. 
It is further transformed into an integral form in section 7 after discussing two relevant technical 
tools, namely, the representations of the Heisenberg-Weyl group (3) and the generalized Fourier 
transform in sections 5 and 6, respectively. An important special case of the Cauchy initial value 
problem for the simple harmonic oscillator is outlined in section 8 and a three parameter general- 
ization of the Fourier transform is introduced in section 9 as a by-product. In sections 10 and 11 we 
solve the initial value problem for the general forced harmonic oscillator in terms of the correspond- 
ing Green function (or Feynman's propagator) and the eigenfunction expansion, respectively, by a 
different method that uses all technical tools developed before in the special case. An extension to 
the case of time-dependent frequency is given in section 12. Then in section 13, we outline impor- 
tant special and limiting cases of the Feynman propagators. Finally in section 14, the motion of a 
charged particle with a spin in uniform magnetic and electric fields that are perpendicular to each 
other is considered as an application; we evaluate a transition amplitude between Landau levels 
under the infiuence of the perpendicular electric field in terms of Charlier polynomials and find 
the corresponding propagator in three dimensions. Solutions of similar diffution-type equations are 
discussed in section 15. 

The Cauchy initial value problem for a forced harmonic oscillator was originally considered by 
Feynman in his path integrals approach to the nonrelativistic quantum mechanics [16], [17], and 
[20] . Since then this problem and its special and limiting cases were discussed by many authors 
fT], [23], [25], [29], [31], [50j the simple harmonic oscillator; [3], [10], [21], [37], [13] the particle in a 
constant external field; see also references therein. It is worth noting that an exact solution of the 
n-dimensional time-dependent Schrodinger equation for certain modified oscillator is found in [30j . 
These simple exactly solvable models may be of interest in a general treatment of the non-linear 
time-dependent Schrodinger equation (see [2S], [21], [36], [S], [IS], [53] and references therein). 
They also provide explicit solutions which can be useful for testing numerical methods of solving 
the time-dependent Schrodinger equation. 




with the Hamiltonian of the form (II. 5p . subject to the initial condition 
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The time-dependent Hamiltonian operator ( 11. 5p has the following structure 

H{t) = Ho + H,{t), (2.1) 

where 

Ho = ^ {aa^ + a^a) (2.2) 

is the Hamiltonian of the harmonic oscillator and 

Hi (t) = H{6{t)a + 6* (t) a^) (2.3) 

is the time-dependent "perturbation" , which corresponds to an external time- dependent force that 
does not depend on the coordinate x (dipole interaction) and a similar velocity-dependent term (see 
[20], [23] and [31] for more details). 

The solution of the stationary Schrodinger equation for the one-dimensional harmonic oscillator 



H,m = E^>, ^o = ^(^-^ + x^J (2.4) 

is a standard textbook problem in quantum mechanics (see [13], [21], [28], [3l], [35], [39], [15], and 
[52] for example). The orthonormal wave functions are given by 

^> = ^>n (x) = -r^== e-^'/^Hr, (x) (2.5) 



with 



oo 



1 Ti ^ TTl 

(X) (X) dx = 5n^ ={ ' ' (2.6) 



where if„ [x] are the Hermite polynomials, a family of the (very) classical orthogonal polynomials 
(see [T], [2], [1], [12], [E], [3S], [Sn], [S], and [IH])- The corresponding oscillator discrete energy 
levels are 

E = E^ = huj\n\^ (n = 0,l,2,... ). (2.7) 



The actions of the creation and annihilation operators (11.31) on the oscillator wave functions (12.51) 
are given by 

a-^n = ^ = Vn + l (2.8) 

These "ladder" equations follow from the differentiation formulas 

-^H^ (x) = 2nHn-i {x) = 2xHn {x) - H^+i (x) , (2.9) 
which are valid for the Hermite polynomials. 
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3. ElGENFUNCTION EXPANSION FOR THE TIME-DEPENDENT SCHRODINGER EQUATION 

In spirit of Dirac's time-dependent perturbation theory in quantum mechanics (see [13], [21], [28], 
[53] , and [13]) we are looking for a solution to the initial value problem in fll.6l) - fll.7l) as an infinite 
series 

oo 

^ = ^{x,t)=Y,Cn {t) ^„ (X) , (3.1) 

n=0 

where (x) are the oscillator wave functions (12. 5p which depend only on the space coordinate x 
and Cn (t) are the yet unknown time- dependent coefficients. Substituting this form of solution into 
the Schrodinger equation (11.61) with the help of the orthogonality property (12.61) and the "ladder" 
relations (12. 8p . we obtain the following linear infinite system 

. dcn (t) f I 



^— ^ = uj\^n+-jcn (t) +6{t)Vn + l c„+i (t) + 6* (t) c„_i (t) (n = 0, 1, 2, ... ) (3.2) 
of the first order ordinary differential equations with c_i (t) = 0. The initial conditions are 

/oo 
^:(x)^o(x) dx (3.3) 
-oo 

due to the initial data (11.71) and the orthogonality property (12. 6p . 

Now we specify the exact form of the function 6 {t) in order to find a particular solution of the 
system (13. 2p in terms of the so-called Charlier polynomials that belong to the classical orthogonal 
polynomials of a discrete variable (see [TT], [12], [H], [38], and [39]). One can easily verify that the 
following Ansatz 

n/2 

C„ (t) = g-i(-(n+l/2)-(n+M))t (e-^^Vn (A)) (3.4) 

vn! 

gives the three term recurrence relation 

Apn (A) = -fiPn+l (A) + (n + /i) pn (A) - upn-i (A) (3.5) 

for the Charlier polynomials Pn (A) = (A) (see [3B] and [22] for example), when we choose 

6{t) = ^ e*{— 1)*, 5* (t) = ^ e-*("-^)* (3.6) 

with the real parameter /i such that < yU < 1. Thus with pn (A) = (A) , equation (13. 4p yields a 
particular solution of the system (13. 2p for any value of the spectral parameter A. 

By the superposition principle the solution to this linear system of ordinary differential equations, 
which satisfies the initial condition (13. 3p . can be constructed as a linear combination 

oo 

Cn (t) = ^ Cnm (t) C„ (0) , (3.7) 

m=0 

where Cnm (t) is a "Green" function, or a particular solution that satisfies the simplest initial con- 
ditions 

Cnm (0) = 6nm- (3.8) 

In the next section we will obtain the function Cnm (t) in terms of the Charlier polynomials; see 
equation (14. 9 p below. In section 5 we establish a relation with the representations of the Heisenberg- 
Weyl group (3) ; see equation (I5.13p . A generalization to an arbitrary function 6 (t) will be given 
later. 
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4. Solution of the Cauchy Problem 

We can now construct the exact solution to the original Cauchy problem in fll.6p - fll.7p for the 
time-dependent Schrodinger equation with the Hamiltonian of the form fl2.1l) - fl2.3p and (13.60 . More 
explicitly, we will solve the following partial differential equation 

1^ = ^ [~'q^ + x ip j + A/2/i I (cos {uj - 1) t) xip + i (sm [u -l)t) — 

subject to the initial condition 

V'(a;,t)li=o = V'o(a;) (-00 < a; < oo) . (4.2) 
By (13.11) . (13.31) . and (13.71) our solution has the form 

00 00 „co 

^ (x, t) = ^ ^„ (x) Cnm (t) / (v) ^0 (v) dy, (4.3) 

n ™_n J —CO 



n=0 m=0 



where 




m °° k 

ml ^-^ k\ 

k=0 

(n+m)/2 

(-I)"-*" g-M ^ g-i(i^(n+l/2)-(n+/.))t 

^/n\m\ 

oo ^ _j4xfc 



X 



fc=0 



in view of the superposition principle and the orthogonality property 

k 



E (k) C (k) e-^^ = ^ (0 < ^ < 1) (4.5) 



k=0 

of the Charlier polynomials (see [2H] for example 



The right-hand side of (14. 4p can be transformed into a single sum with the help of the following 
generating relation for the Charlier polynomials 

oo / \k 

E -^^^ {k) C {k) = e'^^^- (1 - ^^sT (1 - ^Ji2sT (4.6) 

A:=0 

X 2F0 ( -n, -m; ^- 

V (1 - /iiS) (1 - /i2S) 

(See [31], [32], [33], and [22] for more information and [6] for the definition of the generalized 
hypergeometric series.) Choosing /ii = /i2 = and s = e"**//^ we obtain 

C (t) = g-i{(cu-l)n+(n+m)/2)t ^-2^ism^^{t/2) ^-i{ti sin t+{uj/2-fj,)t) y^ 

V2"+"^n!m! 

/ / — \ / 1 \ 
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The hypergeometric series representation for the Charher polynomials is 

<(x)= 2Fo(^-n, -X- --^ (4.8) 

(see for example). Thus 

^ ^-i{tismt+(u>/2~n)t) ^-i{{ui-l)n+{n+m)/2)t i"^-^) 

V2"+'"n!m! 

with P = P (t) = 2A/2/isin (t/2) and, as a result, by substitution of this expression into the series 
( 14. 3p . we obtain the eigenfunction expansion solution to the original Cauchy problem (I4.1l) - fl4.2p . 
We shall be able to find an integral form of this solution in section 7 after discussing representations 
of the Heisenberg-Weyl group N (3) and a generalization of the Fourier transform in the next two 
sections. This complete solution of the particular initial value problem fl4.ip - fl4.2l) will suggest 
a correct form of the Green function (propagator) for the general forced harmonic oscillator in 
sections 10 and 11. 




5. Relation with the Heisenberg-Weyl Group (3) 
Let (3) be the three-dimensional group of the upper triangular real matrices of the form 

= (a,/5,7). (5.1) 

\0 1/ 

The map 

T (a, /?, 7) ^ (x) = e^(^+^^') ^ {x + a) (5.2) 

defines a unitary representation of the Heisenberg-Weyl group N (3) in the space of square integrable 
functions E (— oo, oo) (see [51], [38], and [17] for more details). 

The set {^n(a;)}^o '^^ wave functions of the harmonic oscillator (12. 5p forms a complete 
orthonormal system in {—oo, oo) . The matrix elements of the representation (15. 2p with respect 
to this basis are related to the Gharlier polynomials as follows 

oo 

T («, f3, 7) ^„ (x) = ^ Trnn («, /?, 7) (x) , (5.3) 

where 

POO 

T„,n {a, /5, 7) = / 'i^*m i^) e^^^''""^ {x + a) dx (5.4) 



m 

— oo 



^"^ " -u/2 f ia + f3\"" f icy - [3 



with ly = {a^ + /3^) /2 [38]. A similar integral 



C (^) 



/oo 
i/^ (x + y) [x + z) e-'^'dx = 2"m! ^"""^ L;;"™ (-2yz) (5.5) 
-oo 
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is evaluated in fT5] in terms of the Laguerre polynomials (^) , whose relation with the Charlier 
polynomials is 

c^,(x) = (-/i)-"n!Lr"(/i), (5.6) 
see [38]. Its special case y = z in the form of 

/oo 
(x) Hr. (x) e'^^~y^"dx = 2"m! y""'" L^^" (-2i/') (5.7) 
-oo 

is of a particular interest in this paper. 
The unitary relation 

oo 

X] ("> 7) Tm'n («, /5, T) = ^mm' (5.8) 

n.=0 

holds due to the orthogonality property of the Charlier polynomials (14.51) . 
The relevant special case of these matrix elements is 

•m+n 

T_ (0, /3, 0) = t„„ (/5) = , ^ , , e-^'/^ /2 (n) , (5.9) 

which explicitly acts on the oscillator wave functions as follows 

oo 

e^P-m^ (x) = tmn (/3) *™ (x) . (5.10) 
Relations (gj]), (gS]) and ([53]) imply 

(/32) = ' , , e-(^?-^^^^^^-)/^ (5.11) 

K=0 

X (A + /?2s)"^(/32 + As)" ciin) 
with A = (/?f + /?2 + (s + s^^)) /2, which is an extension of the addition formula 

oo 

J2tmk (A) tnfc (/?2) = tmn (A + /52) (5.12) 
k=0 

for the matrix elements. 

In order to obtain functions c„m {t) in terms of the matrix elements tmn {P) of the representations 
of the Heisenberg-Weyl group, we compare (14.91) and (15. 9p . The result is 

Cnm [t) = (-1)"+"^ e-'^t^sint+{^/2^^)t) g-^{(^-l)n+(n+m)/2)^ ^ (g^^g) 

where P = 2A/2/isin (t/2) . Our solution (14. 3 p takes the form 

oo oo „oo 

(x, t) = E ^„ (x) E c„„ (t) / ^'^(i/)7/^o(l/) (5.14) 

„— n ™_n J —oo 

oo 

^(_l)"e-^*("-i/2)" ^„ (x) 

?i=0 

^ /"OO 

E tmn i/3) / (-l)'" e--*/2 ^m (y) ^0 (y) rfy. 



n=0 m=0 

oo 

-i{tJ, sin t+(Li;/2— 



OO 

X 

TO=0 
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In the next section, we will discuss a generalization of the Fourier transformation, which will allow 
us to transform this multiple series into a single integral form in the section 7; see equations (I7.4p 
and ra- llTTTil below. 

6. The Generalized Fourier Transform 
The Mehler generating function, or the Poisson kernel for Hermite polynomials, is given by 
K, ,) ^ I (.) (,) ^ ^=i= exp (^jyi^^M±fl^ . (6,) 

where '^n {z) are the oscillator wave functions defined by (12. 5p and |r| < 1, r 7^ ±1 (see [H], |41j . 
[19], and [52] for example). Using the orthogonality property (12. 6p one gets 

/oo 
Kr{x,y)^n{y) dy, \r\ < 1. (6.2) 
-00 

Thus the wave functions are also eigenfunctions of an integral operator corresponding to the 
eigenvalues r". 



e-(-/^-)/^ f 2xy-{x' + y')cosT \ 
iCr {x, y) = Ke»r (x, y) = , exp z — (6.3) 



We denote 

V2tt sin T V 2 sin r 

with < r < vr and use the fact that the oscillator wave functions are the eigenfunctions of the 
generalized Fourier transform 

/oo 
ICrix,y)^niy) dy (6.4) 
-00 

corresponding to the eigenvalues e*""^. (See [5], [12] and [H] for more details on the generalized 
Fourier transform, its inversion formula and their extensions. It is worth noting that the classical 
Fourier transform corresponds to the particular value r = 7r/2 [S2]- Its three parameter extension 
will be discussed in section 9.) 

7. An Integral Form of Solution 

Now let us transform the series (I5.14p into a single integral form. With the help of the inversion 
formula for the generalized Fourier transform (see (16.40 with r — r) and the symmetry property 

i/„(-x) = (-l)" Hn{x) 

we get 

/oo 
IC^,/,i-y,z)^miz) dz. (7.1) 
-00 

Then by (15.100 and Fubuni's theorem, 

00 «QQ 

/ (-1)" e""*^' ^™ ^0 (y) dy (7.2) 

POO / POO \ 

/ ( / ^-*/2 (-^' ^) *™ (-2) ^0 (^) 

m=0 -^-^ V>^-oo / 



m=0 

00 
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— oo 
oo 



/C_t/2 (-y, z) ^ tm„ {(3) (z) j ^/'o (y) dyd2 

\m=0 / 

/C-t/2 (-y, ^) (e^^"^„ (2;)) ^0 {y) dydz. 
Now the series (I5.14p takes the form 

CO 

^ [x, t) = e-*('^^i'^*+(W2-A«)t) J2 (-1)" e-^*^'^-^/^)" ^„ (x) 

n=0 

/• /"OO 

X / / /C_t/2 {-y, z) [e'^'^^n {z)) ^0 {y) dydz 




00 

g-*(/.sini+(W2-/.)t) ^ j /C_t/2 (-2/, 2) e*'^" 



X I ^e-^*(--i/2)" M/„ (-x) vI/„ (^) ) ^0 (y) rfyt^^ 

n=0 

g— sin t+{w/2—fi)t) 

f 00 / /"OO 



X / [ / /C_t/2 (-?/, 2;) e^^^/Ct(i/2-^) 2;) rfz J ^0 (y) dy 
J —00 \J —00 J 

in view of the generating relations (16 .ip and (16.30 . Thus 

/oo 
-00 

where we define the kernel as 

/oo 
^t(i/2-u;) (-a;, 2;) e*^^AC„t/2 (-y, 2;) ^2;. 
-00 

This can be evaluated with the help of the familiar elementary integrals 

/•oo \ ~ 

e-^' dx = V^, / e<'^^'+2^^) dz = J- e-^"'!'^ 



(see [9], [I5], and [ID] also). Denoting Ti = t {u — 1/2) and r2 = t/2, from (16.31) we get 
)C_rr {-X, z) e^^^/C_^2 {-y, z) 

cot ri+y^ cot T2)/2 i(l3+x/ sinri+j// sinr2)2 pj(cot ri+cot T2)2^/2 



gj(a;t— tt) /2 



2-7? sin Ti sin r2 

and 



/CX3 
/C_rj (-X, 2;) e*^^/C_T-2 (-2/, -z) c?^ 
-00 



J{ujt-Tr)/2 

cotri+j/ cotT2j/2 



27rv^sin ri sin T2 



X / e 

-00 



00 

i((/9+x/smTi+j// sinT2)2;+(cotri+cotT2)2^/2) 
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As a result 

iM-7r/2)/2 

Qt {x, y) = ™*"^+^ ™*"^)/' (7.9) 

v27rsinu;t 



2i sin cut 



X exp 

with Ti = t (cu — 1/2) , r2 = t/2 and P = P (t) = 2y/2flsm {t/2) . Thus the exphcit form of this 
kernel is given by 

, giM-V2)/2 { (x"^ + y'^)siYiut- ix"^ -y'^)?,iYiiuj -l)t\ ,^ , 

= ■:m^t [ 2.(cos.t-cos(a,-l)i) ) ^'■"'> 

ikf (x, y) 



sin cut (cos cut — cos (tu — 1) t) 



X exp 
where 

kt {x, y) = {x + y) sin {ujt/2) cos ((tu - 1) t/2) - {x - y) cos (tut/2) sin {{cu - 1) t/2) (7.11) 
— A/2/isin (t/2) (cos cut — cos (cu — 1) t) . 
The last expression can be transformed into a somewhat more convenient form 
n, ^ K-w X /sin((cu-l/2)t) sin(t/2) 

^.(x,,)=/C:,(x,,) exp 2. sin cut ^'-''^ 



X exp 



;xsin(t/2) +t/sin((tu - l/2)t)) /? 
i sin cut 



with /? = 2-\/2/Isin (t/2) in terms of the kernel of the generalized Fourier transform (16.31) . Our 
formulas (17.41) and (I7.9p - fl7.12p provide an integral form to the solution of the Cauchy initial value 
problem (14. 1 1) - (14. 2 p in terms of a Green function. 

By choosing t/^q [x) = 6[x — Xq) , where b [x) is the Dirac delta function, we formally obtain 

(x, t) = G (x, xo, t) = e-^('^«i°*+(-/2-A')*)^, (x, xo) , (7.13) 

which is the fundamental solution to the time-dependent Schrodinger equation (14.11) . One can show 
that 

hm (x, t) = V'o (a;) (7.14) 

by methods of [5], [l2] and [52]. The details are left to the reader. 

The time evolution operator for the time-dependent Schrodinger equation (11.61) can formally be 
written as ^ 

[/(t,to) = T fexp [ Hit') dt']] , (7.15) 



where T is the time ordering operator which orders operators with larger times to the left [9], |21j . 
Namely, this unitary operator takes a state at time tg to a state at time t, so that 

^(x,t) = [/(t,to)V^(x,to) (7.16) 

and 

?7(t,to) = t/(t,t')t/(t',to), (7.17) 
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U-'{t,to) = UHt,to)=U{to,t). (7.18) 
We have constructed this time evolution operator exphcitly as the following integral operator 

POO 

U it, to) ^ (x, to) = e-('^^-(*-*«)+(-/2-'^)(*-*°)) / g,^,, (x, y) ^ {y, to) dy (7.19) 



with the kernel given by fl7.9l) - fl7.12p . for the particular form of the time- dependent Hamiltonian in 
f l2.ip - fl2.3l) and (13. 6p . The Green function (propagator) for the general forced harmonic oscillator 
is constructed in section 10 ; see equations fll0.3p - (ll0.8p . 

8. The Cauchy Problem for the Simple Harmonic Oscillator 
In an important special case /i = 0, the initial value problem 



,dip uo f d ip 



i- 



dt 2 V dx"^ 
has the following explicit solution 



+ X if) \ , ip {x,t) l^^Q = ipoi^) {—oo<x< oo) (8.1) 



ip{x., 



OO «QQ 

t) = Yl v[/„ (x) / vl/„ (y) i^y) dy (8.2) 

n J —OO 



n=0 

"OO 



I .{x^ + y'^)cos{ujt) -2xy\ . , . , 
exp I n ■ f ^^ V^o [y) dy. 



^j2Txi sin {yjt) J -oo \ 2 sin {ujt) 

The last relation is valid when < t < tt/u;. Analytic continuation in a larger domain is discussed 
in pj and [50j. 

Equation (18. 2p gives the time evolution operator ( I7.15P for the simple harmonic oscillator in terms 
of the generalized Fourier transform. This result and its extension to a general forced harmonic 
oscillator without the velocity-dependent term in the Hamiltonian are well-known (see [7], [T7], |2Uj . 
|23] . |25j . [33], [SU] and references therein; further generalizations are given in sections 10-12; 
more special cases will be discussed in section 13). 

9. Three Parameter Generalization of the Fourier Transform 

The properties of the time evolution operator in (I7.16p - (l7.19p suggest the following extension of 
the classical Fourier integral 

/■OO 

f{x)= \ Ct{x,y)g{y) dy, (9.1) 



where the kernel given by 



sin {{u - 1/2) t) sin^ (t/2) 2 



A (a;,y) = /C^i (x,?/) exp z — e (9.2) 

\ 2 sm uot J 

/.(xsin (t/2) +?/sin((cu - 1/2) t)) sin (t/2) 

X exp I z E 

\ sin ojt 

depends on the three free parameters t, uj and e. If e = and ut = t we arrive at the kernel of the 
generalized Fourier transform (16.30 . The formal inversion formula is given by 

POO 

g{y)= / C;ix,y)fix) dx. (9.3) 
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The details are left to the reader. Note that, in terms of a distribution, 

C;{x,y)Ct{x,z) dx (9.4) 

Jcot(ujt){y'^-z'^)/2 ^is(z-y) siii{{ui-l/2)t) sui{t /2) / suiLot 
1 



e ' '^^ I' e 

f oo 



2tx sin uot 



e 



ix(z-y)/sinut dx = 5{y-z) 



which gives the corresponding orthogonality property of the £-kernel. These results admit further 
generalizations with the help of the time evolution operators found in sections 10 and 12. 

10. The General Forced Harmonic Oscillator 

Our solution of the initial value problem (14. 1 1) - (14. 2 p obtained in the previous sections admits a 
generalization. The Cauchy problem for the general forced harmonic oscillator 

where / {t) and g (t) are two arbitrary real valued functions of time only (such that the integrals in 
(I10.6p - (ll0.8p below converge and a (0) = 0), with the initial data 

'ip {x,t)\^^Q = ipo (x) (— oo < s < oo) (10.2 

has the following explicit solution 

/oo 

Gix,y,t)^o{y) dy. (10.3 
-oo 

Here the Green function (or Feynman's propagator [T7j, [20], [3l]) is given by 

G (x, y, t) = Go (x, y, t) e*(«(*)-+^'(*)?/+'=W) (10.4 



with 



and 



1 f . {x^ + y )cosujt-2xy \ 

Go [x, y,t)= _ exp U — (10.5 

V2msmujt \ 2 smut J 

1 [' 

a(t) = — / {f (s) sinus + g (s) cos us) ds, (10.6 

sin cut Jq 

h{t)= f^^^i^l^^ds, (10.7 
Jo sincijs 

c{t)= / ( g (s) a (s) - ^a"^ (s)] ds (10.8 







2 



provided a (0) = 6 (0) = c(0) = 0. The case g (t) =0 is discussed in [T7], [20] and [3l| but the 
answers for b (t) and c (t) are given in different forms; we shall elaborate on this later. 



Indeed, the previously found solution fl7.12p - (l7.13p in the special case of the forced oscillator (14.1 
suggests to look for a general Green function in the form (110. 4p . namely, 

^ = M e'^, (10.9 
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where u = Gq {x, y, t) is the fundamental solution of the Schodinger equation for the simple harmonic 
oscillator fl8.ll) and S = a (t) x + b (t) y + c{t) . Its substitution into fllO.ll) gives 

/ da db dc\ ( ^ ^ 2\ \ ^''^ 



where by fllO.Sp 



As a result 



du _ . xcosujt-y ^ (10 11) 
dx sin tut 

da db dc u ^ x cos ut — y /.^.r^N 

—x + —?/ + — = a^ + x/ - -a -{auj-g) — , (10.12) 

dt dt dt 2 smuji 

and equating the coefficients of x, y and 1, we obtain the following system of ordinary differential 
equations 

— {siuLut a (t)) = / (t) smut + g (t) cos tut, (10.13) 

(JjL 

#6(t) = ^^^^ii^, (10.14) 

dt ^ ^ smut ^ ^ 

^c(t) = (7(t)a(t)- V(t), (10.15) 

whose solutions are f ll0.6p - p0.8p . respectively, if the integrals converge. This method is equivalent 
to solving of the quantum mechanical Hamilton-Jacobi equation for the general forced harmonic 
oscillator [3lj. 



Equation (110.71) can be rewritten as 

b{t) = — {sinus a (s)) d cot us — 

Jo Jo 

and integrating by parts 



smcus 



b (t) = — cosut a (t) + / {f (s) cosus — g (s) sinus) ds (10.16) 

Jo 

by (110.131) . With the help of (110. 6p and the addition formulas for trigonometric functions we finally 
arrive at 

b{t) = / {f{s)sinu{s~t)+g{s)cosu{s-t))ds, (10.17) 

sincijt Jq 

which is equivalent to the form obtain in [17], [20] and [31] when g (t) = 0. 



/■* 1 /■* 2 

/ g{s)a{s) ds -\ — / (sinu;sa(s)) d cot us 
Jo 2 Jq 



In a similar fashion, 

c{t) 
and clS 8b result 

1 

c(t) = -sin ut cosut a^ (t) + / sinus a (s) (—f (s) cosus + g (s) sinus) ds. (10.18) 
2 Jo 

This can be transformed into the form given in [T7] and [20] when g (t) = 0. The details are left to 
the reader. 
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Evaluation of elementary integrals results in (17. 121) again in the special case ( 14. ip . The simple 

case f (t) = 2 cos ut and g (t) = gives 

a(t) = ^^^, b{t)=t, c{t) = ^sm2ujt- —t. (10.19) 

The corresponding propagator in (110. 4p does satisfy the Schrodinger equation (110. ip . which can be 
verified by a direct differentiation with the help of a computer algebra system. The details are left 
to the reader. A case of the forced modified oscillator is discussed in l30l. 



11. ElGENFUNCTION EXPANSION FOR THE GENERAL FORCED HARMONIC OSCILLATOR 

Separation of the x and y variables in Feynman's propagator (I10.4I) - (I10.8I) with the help of the 
Mehler generating function (16. ip written as 

oo 

Go (x, y^) = Y^ e-^-(^+i/2)* {x) {y) (11.1) 

fc=0 

gives 

G(x,y,t) = Go(a;,y,t) e^(-+''2^+'=) (11.2) 



k=Q 

oo / oo \ / oo 

e^i'^--'l^)Y^e-'-^' [Y,t^u{a) ^„ (a;) t™^ (6) (y) 



fc=0 \n=0 / \m=0 

oo oo 



= e^^'^--*/^) E E E ^-'^ ^"^'^ 

n=Om=0 \fc=0 

by (IS.lOp . The last series can be summed by using the addition formula (15. lip in the form 



oo 



•m+n 

X (a + 6;2)"(6 + az)" c^'/' (n) , 

with 2 = e~*'^* and = + IP' + 2ab cos cut. As a result we arrive at the following eigenfunction 
expansion of the forced harmonic oscillator propagator 

G (x, y, t) = e'^<^-i-^'-^bsinu^t)/2) g-xV4 (11.4) 

oo oo -n+m 

X E E (^) (^/) ,.n^^ , , + ib + a^)"^ c^^/^ (n) 

n=0 m=0 * 

in terms of the Charlier polynomials. The special case g{t) =0 is discussed in [20] but the 
connection with the Charlier polynomials is not emphasized. 

The solution (110. 3p takes the form 

oo oo „oo 

^ (x, t) = ^ ^„ (x) ^ c„„ (t) / ^^{y)^Q{y) dy, (11.5) 
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where 



Cum (t) = e^(^-("^i°-*)/2) e-^'/^ ^ = (a + (6 + a^)"^ c^/^ (n) (11.6) 

v2"+™-n!m! 

with z = e"*"^* and = + + 2a& cos tut. Functions a = a {t) , b = b {t) and c = c (t) here are 
given by the integrals fll0.6p - fll0.8p . respectively, and lim(^o+ Cnm (t) = ^nm- 

If V'o (a;) = (a^; ''')lt=o ~ (^) ' equation (111.51) becomes 

oo 

'il:{x,t) = 5^c„„(t) ^„(a;). (11.7) 

n=0 

Thus function c„m (t) gives explicitly the quantum mechanical amplitude that the oscillator initially 
in state m is found at time t in state n [SO]. An application to the motion of a charged particle 
with a spin in uniform perpendicular magnetic and electric fields is considered in section 14. 

As a by-product we found the fundamental solution Cnm (t) of the system (13. 2p in terms of the 
Charlier polynomials for an arbitrary complex valued function 6 (t) = (— / (t) + ig (t)) /^/2. The 
explicit solution of the corresponding initial value problem in given by (13. 7p . 

12. Time-Dependent Frequency 

An extension of the Schodinger equation to the case of the forced harmonic oscillator with the 
time-dependent frequency is as follows 

where u (t) > 0, / (t) and g (t) are arbitrary real valued functions of time only. It can be easily 
solved by the substitution 

/■* dr 
T = T{t)= / uj{s) ds, — = uj{t), (12.2) 



Jo dt 
which transforms this equation into a familiar form 

see the original equation ( 110. ip with respect to the new time variable r, where u = 1 and 

/i (r) = 91 (r) = (12.4) 

u [t) OJ [t) 

Therefore by (110.4p -( ri0.5p the propagator has the form 

G (x, y, r) = Go (x, y, r) e^('^(")"+^M?'+^M) (12.5) 

with 

1 /.(x^ + y^) COST — 2a;?/' 



Go {x, y,T) = _ exp i — (12.6) 

V2msmT V 2smr 



and the system (I10.13p -( 110.15p becomes 



^ (sin r a (r)) = /i (r) sin r + gi (r) cos r, (12.7) 
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Thus 



d , , ^ a (t) — Qi (t) , , 

5^=^2_^11^, 12.8 

dr sm r 

±ciT)=g,iT)aiT)-^a'iT). (12.9) 



1 /■* 

a(r) = / {f (s) sin t{s) + g{s) cos t{s)) ds, (12.10) 

SUIT Jq 

■'c.Wa(rW)-jW,^_ (12.11) 



r 



SUIT [S) 

t / 1 



c(r)=/ U(s)a(r(s))--a;(s)a^(r(s)) Ws, (12.12) 



which is an extension of equations fll0.6p - fll0.7p to the case of the forced harmonic oscillator with 
the time-dependent frequency. The solution of the Cauchy initial value problem is given by 

POO 

^{x,t) = / G{x,y,T)^o{y) dy (12.13) 



with r = / u; (s) ds. The details are left to the reader. 



13. Some Special Cases 

The time-dependent Schrodinger equation for the forced harmonic oscillator is usually written in 
the form 

iH—=H^ (13.1) 
dt ^ ^ 

with the following Hamiltonian 

where h is the Planck constant, m is the mass of the particle, u is the classical oscillation frequency, 
F (t) is a uniform in space external force depending on time, function G (t) represents a similar 
velocity-dependent term and p is the linear momentum operator. The initial value problem is 

%n— = - — ^ + -^x^ ^-F{t)x^ + inG (t) — 13.3 
ot 2m ox^ 2 ox 

with 

^ix,t)\,^,^^=^ix,to). (13.4) 

Among important special cases are: the free particle, when uj = F = G = 0; a particle in a 
constant external field, where uj = G = and F = constant; the simple harmonic oscillator 
with F = G = 0. In this section for the benefits of the reader we provide explicit forms for the 
corresponding propagators by taking certain limits in the general solution. 

The usual change of the space variable 
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reduces equation fll3.3p to the form (110.11) with respect to C, with 



/W = ^, 9it) = ^G{t). (13.6) 
The time evolution operator is 

/oo 
G{x,y,t,to) ^{y,to) dy (13.7) 
-oo 

with the propagator of the form 

G (x, y, t, to) = Go (x, y, t, to) e*W*.*o)x+b(t,to)?;+c(i,to))_ (13 g) 

Here 



Go X, y, t, to) = Jtt-t-. T 13.9 

Y zmftsmcu (t — to) 

(imuj 1 1 1 T\ \ 

-— — — {[x +y )cosuj{t- to) - 2xy) 
Insmuj [t — to) 

(^it,to) = — — — - / F s ^ ^ + G s coso; s-to ds, 13.10 

nsmuj [t — to) J to V / 

b{t,to) = -a{to,t) (13.11) 

and 

c{t,to) = (^G{s) a{s,to)-^a^s,to)^ ds. (13.12) 

The simple harmonic oscillator propagator, when F = G = 0, is given by equation (113.91) : see [7], 
[T7] . [20] . [23] , [25] , |29j . [3^ . [50] and references therein for more details. In the limit o; — > we 
obtain 

Go{x,y,t,to) = ,/^^^^^ exp ( 'Z^'^'f, 1 (13.13) 



2mh{t-to) \ 2n{t-to) 

as the free particle propagator [20] . 



For a particle in a constant external field = G = and F = constant. The corresponding 
propagator is given by 

G{x,y,t,to) = ,/^^^ exp ( ^"^(^-^) ) (13.14) 
V ,y, , u; \j2mh{t-to) ^ y2h(t-to) j ^ ' 

f iFix + y) , , iF"^ , 
This case was studied in detail in [3], [10], [20], [21], [37] and [33]. We have corrected a typo in [20] . 
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Figure 1. Magnetic and electric fields in E?. 
14. Motion in Uniform Perpendicular Magnetic and Electric Fields 

14.1. Solution of a Particular Initial Value Problem. A particle with a spin s has also an 
intrinsic magnetic momentum with the operator 



:i4.1) 



where s is the spin operator and is a constant characterizing the particle, which is usually called 
the magnitude of the magnetic momentum. For the motion of a charged particle in uniform magnetic 
H and electric E fields, which are perpendicular to each other (Figure 1), the corresponding three- 
dimensional time-dependent Schrodinger equation 



:i4.2) 



has the Hamiltonian of the form [28] 



H 



1 
2m 



eH 

Px + —y 



s 



:i4.3) 



where p = —ifiV is the linear momentum operator, functions H and F/e are the magnitudes of the 
uniform magnetic and electric fields in z and y directions, respectively. The corresponding vector 
potential A = —yH is defined up to a gauge transformation. Here we follow the original choice 
of [28j (see a remark at the end of this section). 



Since (114. Sp does not contain the other components of the spin, the operator commutes with 
the Hamiltonian H and the 2;-component of the spin is conserved. Thus the operator 's^ can be 
replaced by its eigenvalue = cr in the Hamiltonian 

H = — {v. + —y \ +^pl + ^pl-^H-yF (14.4) 

2m ^ 2m 



1 

2m 



eH 

Px + —y 



with a = —s, —s+1, ... , s— 1, s. Then the spin dependence of the wave function becomes unimportant 
and the wave function in the Schrodinger equation (114.21) can be taken as an ordinary coordinate 



function \E' = \1' 



r,t,a) 
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The Hamiltonian fll4.4p does not contain the coordinates x and z exphcitly. Therefore the 
operators Px and Pz also commute with the Hamiltonian and the x and z components of the linear 
momentum are conserved. The corresponding eigenvalues Px and Pz take all values from — oo to oo; 
see [28] for more details. In this paper we consider the simplest case when the magnetic field H is 
a constant and the electric force F is a function of time t (see Figure 1); a more general case will 
be discussed elsewhere. Then the substitution 

^ (r, t) = e^(^P-+^P--^(*'*o))/^ ^ (y, t) , (14.5) 

where 

S [t, to) = (ll - - + ^ / ^ (^^-^^ 

results in the one-dimensional time-dependent Schrodinger equation of the harmonic oscillator 
driven by an external force in the ^/-direction 

, - F (t) (y - yo) ^ (14.7) 

with 

lei H cpr 

UJH = 2/0 = — ^. 14.8 

mc en 

The Cauchy initial value problem subject to special data 

^ t) 1,^,^ = e^(-f^+^^'^)/^ ^ (y, t,) = e'i^^^+^P^y'' ^{y~ y,) (14.9) 

has the following solution 

/oo 
G{y-yo,yi,tM) ^iv) dv, (14.10) 
-oo 

where the propagator takes the form 

G {y, 7], t, to) = Gi {y, r], t - to) e^W*'*o)?'+^(*'*o)''+^(*'*o)) (14.II) 

with 



Gi{y,V,t) = J—^^^ (14.12) 



( imuju 111 2\ \\ 

^ ^^P TTF^ 7 ((^ + ^ ) coso;i/t - 2yr/) 

\2ftsmco'Ht / 



«(^'^o) = F^ / i"Wsina;H(r-to) rfr, (14.13) 

n Sm LO H [t — iQ) 



to 



b{t,to) = -a{to,t) (14.14) 



and 



c {t, to) = ^ a' (r, to) dr. (14.15) 
See equations (I13.7p - (I13.12I) with G = 0. Function c (t, to) can be written in several different forms. 
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14.2. Landau Levels. In an absence of the external force F = 0, equation (114.71) is formally 

identical to the time-dependent Schrodinger equation for a simple harmonic oscillator with the 
frequency uh- The standard substitution 

^(l/,t) = e-(*-*°)/'^x(l/) (14.16) 
gives the corresponding stationary Schrodinger equation as follows [28] 

2m / 1 \ 

x" + {y - yof] X = 0, (14.17) 

which has the square integrable solutions only when 

e = nujH{n + ]-] , n = 0,1,2,.... (14.18) 



2 



The eigenfunctions are 



X.{y)= \ ^exp(-il^) g„(^). a„ = JJ^, (14.19) 

where if„ {rj) are the Hermite polynomials. 

Thus the total energy levels of a particle in a uniform magnetic field have the form 

En = En{p,,a)=hujH[n+^^+^-^H (n = 0, 1, 2, ... ) . (14.20) 

The first term here gives the discrete energy values corresponding to motion in a plane perpendicular 
to the field. They are called Landau levels. The expression (114.201) does not contain the quantity 
Pxi which takes all real values. Therefore the total energy levels are continuously degenerate. For 
an electron, n/ s = — \e\ H/mc, and formula (114.201) becomes 

En = En (p., ^) = + ^ + ^) + (14-21) 

In this case, there is an additional degeneracy: the levels with n, a = 1/2 and n + 1, a = —1/2 
coincide: En {pz, 1/2) = En+i {pz, -1/2) . 

The three-dimensional wave functions corresponding to the energy levels (114.201) are given by 

vl>„ (r,t,a) = ^n ir,t,p,,pz,a) = e-^"(P--)(*-*o)/^ ^^i^P.+zp.)/n _ 22) 

They are the eigenfunctions of the following set of commuting operators px, Pz, 'Sz, and H with 
F = : 

#^„ = E„^„, s,^„ = a^„, Px^„ = Px^„, p.^n=P.^n. (14.23) 
The orthogonality relation in is 

/ ^*^{r,t,p,,pz,a) mm{r,t,p'^,p'^,a') dxdydz (14.24) 

= {2']Thf 6nm 5aa' S {p^ - p'^) S {pz - p'J , 

where 



27r 

is the Dirac delta function. 



00 



5 («) = ^ / e'°« (14.25) 



00 
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14.3. Transition Amplitudes. In the presence of external force, the quantum mechanical ampli- 
tude of a transition between Landau's levels under the influence of the perpendicular electric field 
can be explicitly found as a special case of our formulas f lll.5l) -( fll.7l) . Indeed, solution (114.101) 
takes the form 

^ (r, t, a) = e-^^(*'*«)/^ ^ ^„ (r, to, cr) ^o) / Xm (v) ^ iv, h) dr, (14.26) 

n=0 m=0 J -oo 

in view of the bilinear generating relation (111.41) . If ip {y.to) = Xm iy) , this equation becomes 

oo 

^ (r, t, a) = e-*^(*'*°)/^ ^ (t, to) (r, to, a) , (14.27) 

n=0 

where coefficients Cnm [t, to) are given by (111.61) in terms of Charlier polynomials as follows 

Cnm (t, to) = e*(^-(^^^(*-*o)-"''^*'^-^^(*~*o))/2) ^-^^4 ^^4_28) 
■n+m 

x^^=== ia + b5nb + a6rct/Hn) 
V2"+'^n!m! 

with 6 = e-''^H(t-to) and -f'^ = a'^ + b'^ + 2ab cos ujh {t - to) . Functions a = a(t,to), 6 = &(t,to) 
and c = c(t,to) are evaluated by the integrals (I14.13l) - (ll4.15p . respectively. The last two formulas 
(I14.27l) - (ll4.28p and (114.61) give us the quantum mechanical amplitude that the particle initially in 
Landau state m is found at time t in state n. For the particle initially in the ground state m = 0, 
the probability to occupy state n at time t is given by the Poisson distribution 

1 

\cnoit,to)f = e-^ ^, fi= -{a^ + b^ + 2abcosuH{t^to)) <l. (14.29) 
nl 2 

The details are left to the reader. 

14.4. Propagator in Three Dimensions. Our particular solutions (114.101) subject to special 
initial data (114. 9p have been constructed above as eigenfunctions of the operators % and pz, whose 
continuous eigenvalues px and Pz vary from — oo to oo. By the superposition principle, one can look 
for a general solution in as a double Fourier integral of the particular solution 



^(r,t)=// aipx^Pz)"^ ir,t,px,Pz) dpxdpz (14.30) 

D 

) 

dpxdpza{px,Pz) e^(-^'^+^^'^)/^ 

D 

/•OO 

xe-''^'^'°^'^ / G{y-yo,r^,tM) ^ iv) dv, 



where functions a {px,Pz) do not depend on time t and S (t, to) is given by (114. 6p . Now we replace 
the special initial data (114. 9p in by the general one 

^{r,t)\,^,^^=<P{x,y,z), (14.31) 

which is independent on p^ (and yo). Letting t — > to in (114. 30p and using the fundamental property 
of the Green function, 

/oo 

G{y-yo,V,t,to) vijl) dr] = <^ {y - yo) , (14.32) 
-oo 
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one gets 

(x, ^) = / / a ip^,p,) ifiy- yo) e^(-f^+^P^)/« dp^dp,, (14.33) 

J J —oo 

where yo is a function of Px in view of (114.81) . Thus 

a (p., P.) V{y- Z/o) = ir^ii^y^ e-*«P-+«^'^)/" didC (14.34) 



(27r^) 

by the inverse of the Fourier transform. Its substitution into (114.301) gives 

^ (r,t) = [ [ dp^dp, e*(-P-+-P-s(*.*o))/?i (14 35) 



oo 

oo 



X / dr] G{y -yo,ri - yo,t,to) 

as a solution of our initial value problem. A familiar integral form of this solution is as follows 

^{r,t)= [ G{r,p,t,to) 0(e,r/,C) d^dr]dC, (14.36) 

where the Green function (propagator) is given as a double Fourier integral 

G (r, p, t, to) = [ [ e'(("-«)P-+("-^)*'^)/^ e-'^m^/" (14.37) 

{2-kH) JJ-oc 

X G {y - yo,ri - yo, t, to) dp^dp;, 

with the help of the Fubini theorem. 

This integral can be evaluated in terms of elementary functions as follows. Integration over p^ 
gives the free particle propagator of a motion in the direction of magnetic field 

Go{z-C,t- to) = exp (^^ (^{z - C) - ll (t - to)) ) dp, (14.38) 



m ( im (z — Q 

exp 



2 



2TTiH{t - to) \ '2fi{t - to) 
by the integral (17. 6p . Thus 

G (r, p, t, to) = exp (^-^ (t - to)) Go{z-C,,t- to) (14.39) 



to 



X G{y - yo,r] -yo,t,to) dp^ 
exp — — (t - to) Go (^ - C, ^ - ^o) 



\ fis 

X Gi (y, r/, t - to) e^(''(*'*«)^+*(*'*o)^+"(*'*°)) 
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exp ( \ X — ^ — : / F (r) dr 



'i{a {t,to) + b{t,to))cp. 



to 



X exp 



eH 



X exp f — ^ (J^ + M + ^) ] tan (uh (t - to) /2) ) dp^ 
\ n \mujH e 

In view of (17.61) . the last integral is given by 

^ exp f ^ ( x - e - ^ / F (r) rfr ) ) (14.40) 



to 



'i{a{t,to) + b{t,to))cp. 



X exp ( — ^ ( + -r7iy + v)Px] tan (a;// (t - to) /2) ) c/px 

/ mujH cot (cu// {t — to) /2) / imwH cot (cu// (t — to) /2) ^2 



where 



with 



e 



= x - e - ^ + r/) tan (cui^ (t - to) /2) + d {t, to) (14.41) 



d {t, to) = -—. 7-—T (14.42) 

eH smojH [t — to) 

X / F (r) (sincu// (r — to) — sinc^H (t — t) — sincj/i' (t — to)) dr. 

J to 

Here we have used (I14.13l) - (ll4.14p . As a result, we arrive at the following factorization of our 
propagator 

G (r, p, t, to) = Goiz-C,t- to) Gi (y, r,, t - to) e'i<t,to)y+mo)n+<t,to)) (14.43) 
X G2 ix,^,y,r],t,to) , 

where Gq {z — C,t — to) is the free particle propagator in (114. 38p . Gi {y, ri,t — to) is the simple har- 
monic oscillator propagator in (114.121) . and 

G2 (x, e, y, V, t, to) = exp (^-^ (t - to)^ (14.44) 

mLUff cot {luh (t — to) /2) / imuou cot {uoh (t — to) /2) 2 

with /3 = /3 (x, e, r/, t, to) given by ffmB -f fm^j ). 

Our propagator can be simplified to a somewhat more convenient form as follows 

G (r, p, t, to) = Go (2; - C, t - to) Gh (x, r/, t - to) (14.45) 

X Gf (x, ^, ?/, 7], t, to) . 
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Here Gq {z, t) is the free particle propagator in the direction of magnetic field. The function 

Gn (x, y, V, t) = exp ^Mnil^tm ^'^'^^^ 

X exp (^-^ (^{{x - 0' + {y- V?) cot {uHt/2) - 2 (x - {y + V^j ) 

is the propagator corresponding to a motion in a plane perpendicular to the magnetic field in an 
absence of the electric field (compare our expression with one in [20], where F = fi = 0, and see a 
remark below in order to establish an identity of two results). The third factor 

G^(x,e,l/,r/,t,to)=e^'^-(*'*°)/^ (14.47) 

with 

Wf (t, to) = n{a [t, to)y + b [t, to)v + c {t, to)) (14.48) 

1 / e 

+ -mujH d {t, to) {d it, to) + 2{x- 0) cot {ujh {t - to) /2) -2—{y + 7]) 
4 \ \e\ 

is a contribution from the electric field. When F = 0, Wp = and Gp = ^■ 

The solution of the Cauchy initial value problem in subject to the general initial data 

^ (r, t)\t^to = ^ ^o) = {x, y, z) (14.49) 

has the form 

^(r,t)= / G'(r,p,t,to) ^(p,to) d^dr^dC,, (14.50) 

which gives explicitly the time evolution operator for a motion of a charged particle in uniform 
perpendicular magnetic and electric fields with a given projection of the spin = a in the direction 
of magnetic field. By choosing ^ (r, to) = 5 {r ~ ro) , where 5 (r) is the Dirac delta function in three 
dimensions, we formally obtain 

^(r,t) = G(r,ro,t,to) (14.51) 

as the wave function at time t of the particle initially located at a point r = Vo- Then equation 
(114.501) gives a general solution by the superposition principle. 

Remark. The vector potential of the uniform magnetic field in the z-direction is defined up to a 
gauge transformation [28] 

A = -yH e,^ A' = A + Vf = -\jH + ]<i:H By = x r (14.52) 
with / {x,y) = xyH/2. The corresponding transformation of the wave function is given by 

VI/ ^ M/' = vl^exp = vl/exp (^^xy)] (14.53) 

and in view of (114.501) 



He J "^V 4/1 



, ^^^^p f ief (x, y) \ ^^^^^ / ief (^, v) \ _ mun (14 54) 

\ he ) \ he J Anih sin {ujHt/ 2) 

X exp (^-^ (^{{x - 0' + {y- V)') cot {uHt/2) - 2^ {xr] - ^y) 



Kl—-x^uj-£i{cosh{2K-l)t)xu + {smh{2K-l)t) — j, (15.1) 
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which is, essentially, equation (3-64) on page 64 of [20j, where we have corrected a typo. The details 
are left to the reader. 

15. Diffusion-Type Equation 

15.1. Special Case. A formal substitution of t —it and ip ^ u into equation (14.11) with u = 2k 
and y/2jl = e yields the following time- dependent diffusion-type equation 

du 
dt 

where the initial condition is 

u (yX,t)\^^Q = Uq (x) (— oo < X < oo) . (15.2) 

As in the case of the time-dependent Schrodinger equation, in order to solve this initial value 
problem, we use the eigenfunction expansion method. Hence the solution is given by 

u (x, t) = ^ ^„ (x) ^ Cnm (t) / (v) uq (y) dy, (15.3) 

n=0 m=0 "^""^ 

where 



X 



by analytic continuation t — >• —it with a; = 2/? and /i = e'^/2 < 1 in (14. 4p and (14.71) . One can easily 
verify that 

lim Cnm [t) = Snm 

and that if < e < a/2 and k > 1/2, k > e^/2, 

lim Cnm{t)=0 (m,n = 0, 1,2, ...) . 

t—*oo 

Thus the limiting distribution is 

lim u (x, t) = {—oo < X < oo) , (15.5) 

t—*oo 

which is independent of the initial data (115.21) . 
Relation (I5.14p becomes 

oo 

, t) = e-(^'/2) -°ht-(K-.V2)t J2 (-1)" e-*(2-i/2)n (15 Q) 



M (X, 

n=0 

X 

m=0 

with /? = /5(t) = -2i£sinh(t/2) . With the help of (Q, (15A0D and the Fubuni theorem we 



OO „QQ 

/ (-1)" e"""*^' ^™ «o (y) 



transform 



J] tn^n (/5) / e-™*/2 vI/„ (_^) (^) rf^ (15.7) 

m=0 





— oo 
oo 
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{-y, Z) I ^ tmn {z) | Uq (?/) c/t/c/z 

\m=0 / 

^e-V2 (-y, 2;) (e'^'^^n (2;)) uo (y) 

where 7 = i/5 = 25 sinh itjl) . The series (115.61) becomes 

00 

n=0 

(e'^^^'n (2)) Mo (y) dydz 

^-(eV2)sinh*-(.-.V2)< /" K ,.,r. {-y , Z) 




X ^ e-*(2- V2)n vj;^ (_^) ^^^^ 

\n=0 / 

_ g-(£2/2) smhi-(K-e2/2)t 

/oo / /"OO \ 
I / (— y, z) e^^i^'g-t(2K-i/2) (—X, z) dz j Mo (y) dy 

-00 \J —00 / 

in view of the generating relation (16. II) . Therefore, the integral form of the solution (I15.3p - fll5.4p is 

/oo 
nt{x,y)uo{y) dy, (15.9) 
-00 

where by the definition 

/oo 
K^-t(2K-i/2) {—X, z) e*^ K^-t/2 {—y, z) dz. (15.10) 
-00 

Denoting ri = q-^C^'^-^/'^) and r2 = e~*/^ we obtain by (16.11) that 

/oo 
Kr,i-x,z)e'''Kr,i-y,z) dz (15.11) 
-00 



7rv/(l-rf)(l-ri) '""^ V 2 (1 - r?) (1 - ri) 

r / (1 - ri^) (1 - rj) 7 - 2ri (1 - rj) x - 2r2 (1 - rf) y 
" i-oo'"^ I (l-r?)(l-ri) " 

(l-r2)(l-r2) ; 

and the integral can be evaluated with the help of an elementary formula 

00 rz" 

'^-'+2''- dz = J- e'''/", a > 0. (15.12) 



(1 + rf) (1 - r2) x^ + (1 - r?) (1 + r|) y^ 



e 

a 
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As a result, an analog of the heat kernel in fll5.9p - fll5.10l) is given by 

(x v) - ^ exp ( - rfrl) jx' + y') + jrj - rl) jx' - y') \ 

^ ' " V-i^-rfr^ V 2 (1 - r?) (1 - r^) ) ^''''^^ 



X exp 



[(1 - rl) (1 - r^) 7 - (ri + (1 - rir2) (x + ?/) - (ri - (1 + rir2) (x 

4(l-rf)(l-ri)(l-r2ri) 



t2 



with ri = e ^'■^'^ ^/^^ r2 = e ^''^ and 7 = 2esinh (t/2) , t > 0. The last expression can be simplified 
to a somewhat more convenient form 

Ht (x, y) = exp (- f (x + y) + ^^—^ {x - y)\ (15.14) 



X exp ( ^ , 2 ^ ) ^nr. (x, y) 



(l -rf)(l- ri) 7' 
1 — r^r 

in terms of the Mehler kernel fl6.ll) . One can show that 



\\m u {x,t) = uq{x) (15.15) 
by methods of [52]. The details are left to the reader. 

A formal substitution of uq (x) = 5 {x — xq) into (115. 9p gives 

U (X, t)=H (X, Xo, t) = e"(^'/2) -'^ht-(.-eV2)t (^5^^g) 

as the fundamental solution of the diffusion equation (115. ip . 
In the limit e — we obtain 

/oo 
K^-2.t {x,y) uoiy) dy (15.17) 
-oo 

as the exact solution of the corresponding initial value problem in terms of the Mehler kernel (16. ip . 
This kernel gives also a familiar expression in statistical mechanics for the density matrix for a 
system consisting of a simple harmonic oscillator [20] . 

15.2. Generalization. A formal substitution of t — > —it and ip ^ u into (110.11) with uj = 2k and 
f ^ f, g ^ —ig yields a diffusion-type equation 

where / {t) and g {t) are real valued functions of time, subject to the initial condition 

M (x, t)|^^Q = Mo (x) (— OO < X < oo) . (15.19) 

The exact solution is 

/oo 
H{x,y,t)uoiy) dy (15.20) 
-oo 

and the Green function can be found in the form 

H (x, y, t) = Ho (x, y, t) e^W^+^'W^'+^W, (15.21) 
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where 

^o(x,y,t) = ^^^^^ -p(^ J (15.22) 

with r = e~2'^*. Indeed, substitution of (115.211) into (115.181) gives the system of equations 

^ (sinh {2Kt) a (t)) = / (t) sinh (2/tt) + ^ (t) cosh {2Kt) , (15.23) 

d , , ^ 2Ka (t) — q it) , ^ 

—b(t) = tVttt^^, 15.24 

dt ^ ' sinh {2Kt) ^ ' 

^c{t) = Ka^ {t) - g{t)a{t) (15.25) 

and the solutions are 

1 /■* 

a (t) = / (/ (s) sinh (2ks) + ^ (s) cosh (2ks)) cis, (15.26) 

Slllll yZiKjL J J Q 

r' 2^(s)-^g(s) ^^^ (15.27) 
Jo smh (2/tsj 

c(t)= [ {na^ (s) - g (s) a (s)) ds (15.28) 

provided a{0) = b (0) = c (0) = 0. 
An analog of the expansion (111.41) is 

oo oo 

H {x, y,t) = Y,Yl ^« ^^/) (1^-29) 

n=0 m=0 

with 

(t) = ^ gC-«t-(afe/2)sinh(2«t)+AV4 /^g gg) 



X (a + 6r)" (6 + ar)"' ai^o 



-m; 



2 
A2 



Here r = e"^'^*^ _\2 _ ^2 _|_ ^2 _|_ 2a5cosh (2Kt) and functions a{t) , b (t) and c (t) are given by the 
integrals (I15.26p - (I15.28I) . respectively. This can be derived by expanding the kernel (115. 2ip in the 
double series in the same fashion as in section 11, or by the substitution t —it, a —ia, b —ib, 
and c — >■ —ic in (lll.4p . The coefficients Cnm (t) are positive when t > 0. 

The solution (115. 20p takes the form 

00 00 „oo 

U{x,t) = Y,'^n{x)J2^nm{t) / m (v) Uq {v) dy . (15.31) 

n=0 m=0 J -oo 

These results can be extended to the case when parameter k is a function of time in equation 
(115.181) . The details are left to the reader. 
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